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A method is proposed to drive an ultrafast non-adiabatic dynamics of an ultracold gas trapped in a 
time- dependent box potential. The resuhing state is free from spurious excitations associated with the 
breakdown of adiabaticity, and preserves the quantum correlations of the initial state up to a scaling factor. 
The process relies on the existence of an adiabatic invariant and the inversion of the dynamical self- similar 
scaling law dictated by it. Its physical implementation generally requires the use of an auxiliary expulsive 
potential. The method is extended to a broad family of interacting many-body systems. As illustrative 
examples we consider the ultrafast expansion of a Tonks-Girardeau gas and of Bose-Einstein condensates in 
different dimensions, where the method exhibits an excellent robustness against different regimes of 
interactions and the features of an experimentally realizable box potential. 



F 



y y ^™ ast good" is a new culinary concept envisioned by the chef F. Adria aimed at creating a diet enjoying of 
two apparently mutually- exclusive features: fast-service and high-quality. When similar ideas are 
invoked in the quantum realm, one faces the adiabatic theorem which imposes a price. The preparation 
of a given target state with high-fidelity, free from spurious excitations, generally demands a long time of 
evolution and the implementation of an adiabatic dynamics. As a result, it comes as no surprise that the recent 
development of shortcuts to adiabaticity (STA)\ has led to a surge of theoreticaP"^° and experimental activity^ 

It is the purpose of this work to show that STA can be implemented in a time- dependent box, the paradigmatic 
model of a quantum piston The relevance of this confinement is enhanced by the development of experimental 
techniques to create optical billiards for ultracold gases^^'^^, paint arbitrary potential shapes for Bose-Einstein 
condensates^^, and realize all-optical boxes^^ and analogous traps in atom chips^^. It has the potential to greatly 
advance the field of quantum simulation, facilitating the connection between ultracold atom experiments and 
condensed matter systems. 

The prospects of finding STA for box confinement seem challenging at the very least, based on the following 
considerations. Since the early insight by M. Moshinsky in 1952, it is known that the rapid expansion of matter- 
waves initially localized in a region of space exhibits quantum transients and excitations in the form of density 
ripples, a phenomenon referred to as diffraction in time^°"^^. Moreover, reflections from the walls of the trap 
generally lead to Talbot oscillations and the formation of a quantum carpet in the time evolution of the density 
profile^^'^'*. On top of that, a series of work during the last decades have shown that the suppression of excitations 
in an expanding box is inevitably constrained by the adiabatic theorem both in the non-interacting^^"^^ and mean- 
field regime^^. 

In spite of these results preventing ultrafast excitation-free dynamics in box traps, we show that fast non- 
adiabatic expansion or compression is in fact possible, allowing preparation in a predetermined finite time of the 
same target state than the adiabatic evolution. The key to this shortcut to adiabaticity is the use of an auxiliary 
potential term superimposed on the box trap. 

We shall start by introducing a dynamical invariant in a time -dependent box trap at the single-particle level, 
and use it to derive a scaling law that governs the time evolution. The inversion of this scaling law will be the key to 
engineer the STA. This result will be extended to many-body systems with a broad family of interactions. Finally, 
we shall illustrate the method with some examples relevant to experiments with ultracold gases, showing that the 
STA is robust in the presence of interactions and experimental imperfections. 

Results 

Dynamical invariants and self-similar dynamics. For a time- dependent box of width ^(t), the scaling laws 
governing the dynamics of the expanding eigenstates reported to date are associated with trajectories of the form 
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Figure 1 | Time-modulation of the potential trap along a shortcut to 
adiabaticity. (A) Evolution of the width ^(t) of a time-dependent box 
during a shortcut to different adiabatic expansions and a compression in a 
finite time t. (B) Frequency of the auxiliary external harmonic potential 
required to assist the self-similar dynamics, changing character from 
expulsive to attractive along the expansion process. The opposite sequence 
is required for a compression. The transient expulsive potential is 
responsible for the speed-up of the process. 

^{t) = [af + ht + cf (with a, h, c real constants) which turn out 
to be unsuited for engineering a STA (see below). Nonetheless, given 
a time- dependent Hamiltonian H(t), it is possible to build a 
dynamical invariant^^, X{t) such that 



dX{t) dl{t) 1 



dt 



dt 



+ -[i{t)Mt)]=o. 

m 



with spectral decomposition i(t)= ^n\4'n{'^)){^n{'^)\ terms 
of the set of eigenmodes |(/>n(0) with eigenvalues An- This is a 
particularly useful basis to describe the time evolution of an initial 
state ^, by the superposition ^{x,t) = 5I]nexp(/a„)|(/>„(0), where 
the Lewis -Riesenfeld phase is given by a„ = {(l)„{t')\ihdt' —H{t')\ 
(l)^{t'))dt' /h, and can be understood as the sum of the dynamical 
phase and the Aharanov-Anandan phase^^. For a time- dependent 
box of width ^(t) and initial width ^(0) = ^q, a dynamical in- 
variant exists^ ^ 



1 = 



1 e{t) 

2m a 



m 

-m-—-x 



with eigenvectors {x\(l)^{t)) = [2/^{t)Yexp 



and eigenvalues In = — -y with kn 
2m 



(2) 



sm[nnx/^{t) 



nn/^Q. The Lewis -Riesenfeld 



phase can be computed to be (Xn{t)-- 



rj{t) with 



f]{t)= ^Qdt'^l/^^{t'). The condition in Eq. (1) for X to be an 
invariant requires the box potential to be supplemented with an 
auxiliary harmonic term 



U^'^ix.t) -- 



1 



m 2 



(3) 



whose frequency Q{t) = 



m 
m 



is dictated by the ratio of the 



acceleration of the trajectory ^{t)^ and the trajectory itself. Note 
that ^{t) > 0, so that for ^{t)>0, Q{t) is purely imaginary, and 
the auxiliary term L^'''''' is a repulsive harmonic potential. If t{0) = 0, 



[H(0),X(0)] =0 so H(0) andX(O) have common eigenstates. Fur- 
ther, if '^(0)=0 holds as well, then U''''''{x, 0) = 0 and an eigenstate 
^„(x, 0) of the box att = 0 evolves into ^„(x, t) = exp[/a„(0](/>n(Oj 
a key observation to engineer a STA as we shall see. We note that 
the experimental implementation of V'^ix, t) can be assisted by the 
same techniques used to create the box potential: the use of a blue- 
detuned laser^^ or direct painting of the required trap^^"^^. 

Shortcuts to adiabaticity: inverting the dynamical scaUng law. We 

next discuss how to implement a non- adiabatic expansion of the 
box by a factor y(T) = ^(t)/(^o in a given finite-time t suppressing 
excitations in the final state. We shall impose the condition U''''''(x, 0) 
= U^'^ix.T) = 0. As in the adiabatic case, in a STA the time evolution 
of an eigenmode of the initial box should reproduce an eigenmode of 
the final trap. As at t = 0, this can be enforced by imposing the 
condition c^(t) = (J(t) = 0. The set of boundary con- 
ditions at f = 0, T excludes the possibility of a linear ramp, as well 
as the family of trajectories, ^{t) = [af -\-bt-\- cf, considered so far in 
the literature^^"^^ However, it suffices to determine a polynomial 
ansatz for the trajectory ^(t) = a^t^, i.e. a scaling factor of the 

form y(^)= ^ = l + [y(T)-l]s^[10 + 3s(2s-5)], with s = tlx. 
so 

This further determines the required time-dependent frequency of 
the auxiliary harmonic potential U^'^iXy t) according to Eq. (3), 



y(T) 



2y(t) 



-60s[l+s(2s-3) 



(4) 



The trajectory, displayed in Fig. 1 shows that during an expansion 
U^'^iXyt) becomes an expulsive potential in an early stage (t < t/2), 
providing the speed-up required to achieve the STA in an arbitrary 
finite time t (bounds in the presence of perturbations will be discussed 
below). In a subsequent stage, t > t/2, Q\t) changes sign and U^'^(x, t) 
becomes a trapping potential, slowing down the expanding mode and 
reducing it to an eigenstate of the final Hamiltonian at t = t. Pre- 
cisely the opposite behavior is exhibited during a fast nonadiabatic 
compression. Provided that an arbitrary Q^(t) dependence can be 
implemented, a STA has no lower bound for t (notice however that 
Q(t) ~ T"^). By contrast, the adiabatic condition 



max max 

t n,k 



{E„{t)-Ek{t)) 



<1, 



(5) 



leads to the requirement mnk^{t)^{t) /h^l- 
Note that the energy of the expanding mode 



{H{t))^=En{t) + 



mi jt) 
12 



3 



(6) 



has two contributions, the first one being the adiabatic energy E^it) 
= En(0)i(Oy/i(ty, and the second one depending explicitly on c^, so 
that (H(t))^ =En{T) given that a STA demands c^(t) = 0. This rela- 
tion illustrates the fact that a STA is associated with a non -adiabatic 
evolution, which reproduces the adiabatic result at the end of the 
process. Moreover, STA work as well for excited states: the time 
evolution of the n-th eigenstate of the initial trap leads to the n-th 
eigenstate of the final trap at t = t. As a result, STA in boxes pave the 
way for fast population-preserving cooling in the following sense. 
Given a system described by the canonical ensemble with a density 

matrix e~^^y/Tr , where P = l/Zc^T, ks is the Boltzmann 

constant and T is the temperature, the final temperature reads 



r(T) 



m 



(7) 



These results strictly hold for a box with infinite walls at x e [0, 
i{t)]. However, in the all-optical trap reported in^^, the end-cap lasers 
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Figure 2 | Space-time density plot of the scaled density profile n{xyt)y{ t) during a shortcut to an adiabatic expansion for the ground- state mode of a box 
with smooth boundaries. A high fidelity-state results from the STA (left) in a box created with Gaussian barriers of height Vq = 10^^^/ m^l and a = ^o/20 
(y = 10). In the absence of W'^'^ix) (center), and for a linear ramp box width ^(t) (right) the cloud is excited and exhibits fringes associated with the 
reflections of the boundary. The color coding varies from purple to red as the function value increases. The dashed line follows the trajectory of the moving 
walU(t). 



providing the box walls have a Gaussian profile which smooths the 
potential in a length scale cr, i.e. a box trap of the form 



exp 



X 

2^ 



+ exp 



2<t2 



(8) 



A similar smoothing occurs in other physical realizations^^. Since this 
smoothing is expected to be the most significant deviation of labor- 
atory potentials from the ideal infinite box, it is important to consider 
its effect on the self-similar dynamics required for a STA. This can be 
quantified by the overlap between the states resulting from the 
expansion in an idealized and realistic box trap, and 
\^na(t)) respectively. Glearly the role of cr decreases (increases) dur- 
ing a expansion (contraction) of the box. The numerical solution of 
the time- dependent Schrodinger equation for this box potential is 
shown in Fig. 2, where the STA is compared with both the poly- 
nomial and linear expansion of the box in the absence of U^''''(x,t). 
It is seen that the STA is the only successful strategy and that the 
process is robust even for a substantial smoothing of the potential 
barriers, where the target states deviate from those of an idealized 
box. 

We have further explored numerically the dynamics under "con- 
catenated STA", in which the overall expansion is splitted into a 
sequence of k STA with either constant expansion factor y or con- 
stant box size increment between consecutive steps. The efficiency of 
the process exhibits a non-monotonic improvement with increasing 
k, suggesting a natural scenario where STA techniques could be 
combined with optimal quantum controP'^. 

Beyond adiabatic invariants: Shortcuts to adiabaticity in inte- 
racting many-body systems. Knowledge of the adiabatic invariants 
for a single particle in a time- dependent box has provided us with the 
insight of assisting the dynamics with an auxiliary potential to design 
STA for expansions and compressions in a finite time, without in- 
ducing quantum transients associated with diffraction in time, which 
are ubiquitous in this type of scenario^^. The technique can be directly 
applied to non- interacting gases and other many-body quantum fluids 
which can be mapped to non-interacting systems. It further suppresses 
the Talbot dynamics associated with quantum carpets woven by the 
density profile typically observed in boxes^^ and the question naturally 
arises as to its applicability to interacting systems^^. The presence of 
interactions, e.g. a two-body potential, hinders the exploitation of the 
superposition principle in terms of the eigenmodes of the Lewis- 
Riesenfeld invariant. However, we note that to design a STA it suf- 
fices to enforce a self- similar dynamics and ultimately no knowledge 
of adiabatic invariants is required. As a result, we next consider a 
broader family of many-body systems, confined in a box, defined by 
the Hamiltonian 



2m 



Aq,+ t/»-(q,,() 



+eE^(q,-q;) (9) 



where q, e IR , Aq, is the Laplace operator in dimension D, the 
auxiliary term is now given by 



U'"^{q,t) = 



1 ^(0, ,2 



(10) 



and the two-body interaction potential obeys V(>^q) = yl"°'V(q), e.g. for 
the Fermi-Huang pseudopotential describing s-wave scattering in 
ultracold gases, a = D. For a hard- wall box, = |q|/ e [0,(^(0]; we 
shall relax below this approximation and consider realistic potential 
boxes as those created in all- optical setups. The case D = 1 
corresponds to a box with one-wall moving (the symmetric case in 
which both walls move in opposite directions can be obtained by a 
Duru transformation^^). For D = 2,3 cyhndrical and spherical 
symmetry is assumed respectively. Without loss of generality, we 
choose the dimensionless time- dependent coupling constant e = e{t) 
to satisfy e(0) = 1. A stationary state ^{t) = ^(qi, q^;0 of N 
particles and chemical potential /n follows for ^ > 0 the evolution 



^(t) = y-— exp 



2yh 



qi 

y{t) 



,(11) 



with the boundary conditions ^{t) = 0 for |q|/ = ^{t) {i = l,...,N, in 
addition to ^(t) = 0 for |q|/ = 0 in D = 1), as long as 



<t)=y{t) 



(12) 



which can be implemented exploiting a Feshbach resonance or 
modulating the transverse confinement in anisotropic systems'*'^^. 

The self- similar dynamics in a STA leads to a scaling of all local 
correlation functions. In particular the density of a given many-body 
state follows the law n{q,t) = j . . . <^q]vl^(q'q2' • • • '^^5 01^ ^ 
n[q/y{t),0]/y{t)^ . By contrast, non-local correlation functions 
exhibit a non-trivial dynamics. The one-body reduced density matrix 
^i(q,q'; t) =Njdq2 . . . dq^^iq.q^ . . . ,q^; 0^*(q',q2 • • • An'^ t) of 
a state obeying Eq. (11), follows the scaling law gi{q,q^;t) = 



exp 



2yn 



^i(q/y,qyy; 0), analogous to that observed 



under harmonic confinement^^ The additional phase factor induces 
a major distortion of the momentum distribution, n{k,t) = 
J(iq(iq'^i(q,q'; t) exp [/k-(q — q')]/ (2n)^. However, a STA ensures 
that dXt = T the Lewis -Riesenfeld phase factor vanishes, so that the 
final state exhibits the same correlations of the initial state scaled by a 
factor y(T), 
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^i(q,q';T) = 

y 

n(k,T) =y(T)^n(y(T)k,0) 



;0 



(13) 



Examples. In the following we shall illustrate different aspects of 
shortcuts to adiabaticity in some paradigmatic models. 

We shall first consider the evolution of correlations in a one- 
dimensional cloud of ultracold bosons in the limit of hard-core con- 
tact interactions, this is, in the Tonks- Girardeau (TG) regime^^. 
This system, as well as its lattice-version, has become a favorite 
test-bed to study the breakdown of thermalization and adiabaticity^^. 
Its many-body ground state is given by the Bose-Fermi mapping^^, 
1 



^Tciqu-'-qN)-- 



l<j<k< 



N^{qk 



-qj) det 

j,k— 1 



where 



e{q) = 1 (-1) if jc > 0 (< 0) and e(0) = 0. In a STA, the self-similar 
dynamics is inherited from the single-particle orbitals ^j(qh t) 
whence it follows that no tuning of interactions is required. Its den- 
sity exhibits a scaling law for all t. The same holds true for the 
entanglement entropy with respect to a bipartition [0, a^(t)] (a < 
ly^. The self-similar dynamics breaks down for the momentum 
distribution, which can be computed efficiently^^, and we shall focus 
on its evolution along a STA. Different snapshots are depicted in 
Fig. 3A, and confirm that during an expansion the cloud is acceler- 
ated during the interval [0, t/2] and slowed down during [t/2, t]. The 
reverse sequence, is observed in a fast frictionless compression. The 
axis are scaled up by the expansion factor y(0 in such a way that for 
an adiabatic dynamics, curves at different times would collapse into a 
single curve. Along a STA, the width and mean of the momentum 
distribution do not remain constant and change along the process. 

A similar distortion of correlations, known as dynamical fer- 
mionization, occurs in the dynamics of a cloud suddenly released 
from an arbitrary trap^^. Under ballistic dynamics the asymptotic 
momentum distribution in a ID expansion evolves to that of the 
dual system, a spin polarized Fermi gas. In particular for a cloud 
released from a box the exact time evolution is not self-similar'*^ 
but dynamical fermionization is observed^*'^^. However, under a 
self-similar scaling law, the asymptotic n{k) maps to the density 
profile of the initial state^^ and no dynamical fermionization 
occurs. This is the case of relevance to STA, where the dynamical 
scaling law in Eq. (11) holds. (We note that the case of the initial 
harmonic confinement is singular in that the free expansion is self- 
similar and that the the single-particle eigenstates can be written in 
terms of Hermite polynomials, which are eigenfunctions of the 
continuous Fourier transform. As a result the asymptotic 
momentum distribution can be related to both the initial density 
profile and the momentum distrubtion of non-interacting fer- 
mions. See [5] for a discussion of STA in harmonic traps.) 
Moreover, this distortion of correlations is not restricted to expan- 
sion processes. Along a STA, this is shown in Fig. 3 for both 
expansions (A) and compressions (B). This is a spurious effect 
for the purpose of STA, which is to reproduce the adiabatic result 
in a finite short time. Indeed, the distortion induced during the 
first half of the STA associated with the accelerated expansion or 
compression, is compensated in the second half of the dynamics, 
in such a way that the correlations of the initial state are recon- 
structed at t = T and scaled by a factor y(T). 

We next turn our attention to the design of STA for a BEC in time- 
dependent box trap, where different strategies can be adopted 
depending on the dimensionality and the regime of interactions. 
The time-dependent Gross-Pitaevski equation (TDGPE) governs 
the evolution of the normalized condensate wavefunction 0(q, t), 
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Figure 3 | Time evolution of the momentum distribution of a Tonks- 
Girardeau gas confined in box during an ultrafast expansion and 
compression providing a shortcut to adiabaticity. In a shortcut to an 
adiabatic expansion (A), due to the auxiliary external potential in Eq. (3), 
the cloud is accelerated in the first half of the process. At s = 1/2 there is a 
partial mapping of the momentum distribution to the density profile of the 
initial state of the cloud. In the second half of the process, the cloud is 
slowed down and at t = t it matches the result from an adiabatic dynamics: 
the state exhibits exactly the same correlations of the initial state, scaled up 
by the expansion factor 7(1) = 10 (N = 10). The reverse sequence is 
observed in a shortcut to an adiabatic compression (B), where the target 
state is reached with high-fidelity already for s = t/z = 0.75, after which the 
dynamics is actually adiabatic (7(1) = 1/10). 



(14) 



ihdM<l, t)=[- £A,+ t/--(q, t)+gomq, t)\']<P{q, t), 

|q| e [0,^(0], 

for which adiabaticity conditions have been reported^^. The ansatz 



(^{t)=y 2 exp 



■.mjqjj .fir]{t) 
I 1 



2yn h 
satisfies the TDGPE provided that 



(I)[q/y(f); 0], (15) 



yl-D ■ 



(16) 



These relations constitute the box analogue of the well-known 
Castin-Dum-Kagan-Surkov-Shlyapnikov relations in harmonic 
traps^^'^^. The two-dimensional case is special since the scaling law 
holds when ^2d(0 is kept constant. 
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Figure 4 | Shortcut to an adiabatic expansion of a Bose-Einstein condensate. The upper row shows the space-time evolution of a quasi- ID BEC 
(n(xyt)y(t)) while the lower row corresponds to a quasi-2D {n{r,t)y{tY) BEC cloud. The shortcut to adiabaticity (left) is compared with a polynomial 
(center) and a linear ramp (right) of the box width ^(t) in the absence of the auxiliary harmonic potential (gI^q = 1/20, = 10/z^/ m^l~^, and y = 10). 



Figure 4 is a set of numerical solutions of the time- dependent 
Gross-Pitaeveskii equation that illustrate the robustness of the STA 
for realistic BEC experiments. We consider a box trap with Gaussian 
barriers and in all numerical simulations interactions are kept con- 
stant, i.e. gM = goiO), deviating from the ideal prescription in Eq. 
(16). The top row illustrates the dynamics for a quasi- ID BEC. The 
(one-body) fidelity between the resulting state 0(t) and the ground 

state of the final box is J^y ^ lo = | (^^{t) \^gs^\ = 0.91 1. For smal- 
ler values of y, the fidelity is even higher Ty=^ = 0.999 as expected, 
given that implementation of the exact STA requires a smaller tuning 
of ^1. The bottom row shows the dynamics of a quasi-2D cloud, 
which requires no interaction tuning in a STA, but is more sensitive 
to the smoothness of the box boundaries, Ty^ lo = 0.988. 

It is noteworthy that in the Thomas -Fermi regime, the kinetic term 
contribution can be neglected and it is possible to induce an exact 
self-similar dynamics (and a STA) exclusively with the help of an 
external field. Then, the scaling ansatz is a solution of the TDGPE 
provided that 

^W=f^, ^\t)=-\. gD{t)=go{0). (17) 
J y(0 C 

This regime is particularly robust against the smooth boundaries of 
physically realizable box potentials. The simulations correspond to 
the most delicate regime with moderate mean-field interactions, both 
far from the non-interacting and Thomas-Fermi limits. 

Discussion 

In conclusion, we have presented a method to drive an ultrafast 
dynamics in a time- dependent box trap which reproduces the adia- 
batic result at the end of the evolution. The method is assisted by an 
auxiliary external harmonic potential which provides the speed-up 
and is applicable to a large family of both non-interacting and inter- 
acting many-body systems supporting dynamical scaling laws, where 
it not only leads to a robust expansion of the density but also pre- 
serves the non-local correlation functions of the initial state, up to an 
expansion factor. The proposal is applicable to realistic box poten- 
tials and can be implemented in the laboratory with well-established 
technology. Its applications range over all scenarios requiring a 



shortcut to adiabaticity, i.e., probing strongly correlated phases, pre- 
venting decoherence, the effect of perturbations and atomic losses. 
The method can be directly applied as well to ultrafast population- 
preserving cooling methods, quantum heat engines and refrigera- 
tors'*^ providing an alternative to the paradigmatic model of a 
quantum piston*^. 
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